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Abstract. We present our current studies and our future plans on microscopic potential based on effective 
nucleon-nucleon interaction and many-body theory. This framework treats in an unified way nuclear struc¬ 
ture and reaction. It offers the opportunity to link the underlying effective interaction to nucleon scattering 
observables. The more consistently connected to a variety of reaction and structure experimental data the 
framework will be, the more constrained effective interaction will be. As a proof of concept, we present some 
recent results for both neutron and proton scattered from spherical target nucleus, namely '*°Ca, using 
the Gogny DIS interaction. Possible fruitful crosstalks between microscopic potential, phenomenological 
potential and effective interaction are exposed. We then draw some prospective plans for the forthcoming 
years including scattering from spherical nuclei experiencing pairing correlations, scattering from axially 
deformed nuclei, and new effective interaction with reaction constraints. 
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1 Introduction 

Producing satisfactory data evaluations based solely on 
many-body theories and effective nucleon-nucleon (NN) 
interaction is a long term project. Two keys to success are: 
(i) robust and well tested nuclear reaction codes, such as 
TALYS [T] or EMPIRE [5], flexible enough to incorpo¬ 
rate new microscopic models and (ii) microscopic inputs 
such as optical model potentials, nuclear level densities, 
7 -ray strength functions, and fission properties, all based 
on effective NN interaction. Advances made along this line 
provide more opportunities to connect the effective NN in¬ 
teraction to a broad body of structure and reaction data 
and as a matter of fact to improve its parametrization. 


deformations determined within axially deformed HFB 
with Gogny DIS interaction all over the nuclide chart [S]. 
Effective theories are usually based on phenomenological 
parametrizations of the effective NN interaction, such as 
Skyrme [51[7] or Gogny forces [SII^ ITIJIITT] . In the follow¬ 
ing we mainly focus on Gogny interaction. The fit of the 
interaction mostly relies on connection made with struc¬ 
ture data through the effective theory. Those constraints 
are then completed by physical filters coming from infi¬ 
nite matter calculations. Up to now, two strategies have 
been adopted. The first and original one uses a restricted 
HF model where single-particle orbitals are approximated 
by harmonic-oscillator wave functions for simplicity. This 
makes possible to determine parameter sets of the inter¬ 
action from a limited number of constraints by matrix 
inversion and from physical filters. This strategy has been 
applied with success to the determination of D1 [3], DIS 
[lO] and DIN m versions of Gogny interaction. Recent 
multiparticle-multihole configuration mixing studies m 
na have motivated the elaboration of a generalized Gogny 
ties and (quasiparticle-)random-phase approximation ((Q)RliAl)eraction with finite-range density, spin-orbit and tensor 


1.1 Nuclear energy density functional 

Many-body theories based on effective NN interaction such 
as Hartree-Fock(-Bogolyubov) (HF(B)) for static proper- 


[3] or five-dimension collective Hamiltonian (5DGH) 
for dynamical properties have proven their ability to de¬ 
scribe a wide range of nuclear structure observables, in¬ 
cluding binding energy, charge radius, deformation, exci¬ 
tation spectrum, density and spectroscopic factors, this for 
nuclear masses with A > 5. As an illustration of the ex¬ 
tended reach of the method, we show in Fig. [1] nuclear 
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terms. Along that line, Gogny D2 interaction with a finite 
range density term has been designed M- The second 
strategy is more based on brute force using HFB with a 
self-consistent account of quadrupole correlations energies 
within the 5DCH approach. This strategy has been used 
to develop the DIM parametrization of Gogny interaction 
which reaches a rms deviation with respect to the 2149 
measured masses of only 798 keV [8]. This improvement 
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Fig. 1. Ghart of the nuclides showing deformations obtained with axially deformed HFB/DIS [5]. 


of the parametrization has been done conserving as much 
as possible the virtues of former Gogny interactions. 

As a consequence, the global character of effective theories 
as well as their accuracy and their relatively low compu¬ 
tational cost make them well suited to fulfill the needs for 
accurate nuclear data files in a reasonable time scale. 


1.2 Nuclear codes ingredients based on effective 
interaction 

Various phenomenological ingredients of reaction codes 
have been replaced by their microscopic counter parts 
based on effective interaction. Connecting the dots, it will 
be possible to obtain satisfactory evaluations on the ba¬ 
sis of reliable and accurate microscopic inputs only m- 
A good example of this on going work is the use made 
of results on ^38^ isotope described within QRPA with 
Gogny force [Min]- Those results have made possible a 
microscopic description of preequilibrium without ad-hoc 
“pseudo-state” prescription [18]; see also Dupuis et al. in 
this issue. Along the same line, nuclear level densities have 
been obtained using a temperature-dependent HFB ap¬ 
proach with Gogny interaction US]. 7 —ray strength func¬ 
tion studies based on QRPA and Gogny interaction have 
also been developed [201I2I]. Concerning optical potentials, 
TALYS uses phenomenological potentials mm in the 
regions where data are available and a modified version of 
the semi-microscopic Jeukenne-Lejeune-Mahaux |24] else¬ 
where. Neither of those methods allows a direct connec¬ 
tion with NN interaction. In order to fulfill this lack, we 
are interested in developing a microscopic potential based 
on the effective NN interaction. 


1.3 Microscopic and ab-inito potentials 

Depending on the projectile energy and the target mass, 
various strategies have been adopted in order to deal with 
elastic scattering starting from NN interaction. We now 
expose the pros and cons of those different methods. In 


the following, ab-initio (microscopic) refers to methods 
based on bare (effective) NN interaction. 

Nuclear matter models [25] provide reasonable descrip¬ 
tions of nucleon elastic scattering at incident energies above 
about 50 MeV [26], even up to ^1 GeV (27]. The method 
is based on the folding of a matter density and a g-matrix 
effective interaction built from bare NN interaction. The 
density is provided by an ad-hoc prescription such as mean- 
field or beyond mean-field approaches. Work toward a 
consistent treatment of both density and g-matrix is in 
progress [55]. Recent ab-initio calculations address the 
issue of reactions involving light nuclei and low-energy 
regime. The resonating group method within the no-core 
shell model, has successfully described nucleon scattering 
from light nuclei [29] . This ab-initio model has recently 
been extended to include three-nucleon forces for nucleon 
scattering from "‘He mm- They deal with '^He, "‘He and 
‘“Be targets and incident energies below 15 MeV. An¬ 
other method, the Green’s function Monte Carlo method 
has been used to describe nucleon scattering from "‘He 
in particular the phaseshift of the first partial waves be¬ 
low 5 MeV incident energy ^5]. Other ab-initio calcula¬ 
tions handle magic nuclei. Among them, the self-consistent 
Green’s function (SCGF) method has been applied to op¬ 
tical potential calculations for 40.48.60(7g^ targets [33]. The 
SCGF potential is compared with a phenomenological dis¬ 
persive potential [34]. This model underestimates nuclear 
radii and, as a consequence, is not well suited for scat¬ 
tering calculations. Further studies including three-body 
forces may cure this issue. Moreover work on Gorkov- 
Green’s function theory is in progress to extend SCGF 
to nuclei around closed-shell nuclei [331133] . Finally, the 
coupled-cluster theory has been applied to proton elastic 
scattering from "‘“Ca [37]. Cross section at 9.6 MeV and 
12.44 MeV center-of-mass energy are compared with data. 
They observe a lack of absorption. 

Although ab-initio methods have made progresses in han¬ 
dling light and magic-nuclei, they are still yet suited nei¬ 
ther for heavy targets nor for high incident energy projec¬ 
tiles. Another option is to build the potential starting from 
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an effective NN interaction. The price to pay is to break 
the explicit link with bare NN interaction. The advantage 
is once again the extended reach of effective theories and 
the wealth of results already available. 

The so-called nuclear structure method (NSM) for scatter¬ 
ing [551[5^HniHT] relies on the self-consistent HF and RPA 
approximations to the microscopic optical potential |42] . 
The former is a mean-field potential; the latter is a po¬ 
larization potential built from target nucleus excitations. 
This method applies to double-closed shell spherical tar¬ 
get nuclei well described with RPA. A strictly equivalent 
method, the continuum particle-vibration coupling using 
a Skyrme interaction, has been recently applied to neutron 
scattering from [13]. They neglect part of the residual 
interaction in the coupling vertices. In addition, they do 
not address the issue of the double counting of the uncor¬ 
related second-order diagram. Other approaches aiming 
at fitting a Skyrme effective interactions including reac¬ 
tion constraints are in progress, where optical potential is 
approximated as the HF term and the imaginary part of 
the uncorrelated particle-hole potential neglecting collec¬ 
tivity of target excited states miis]. A recent application 
of NSM with Gogny interaction is presented in Ref. m- 
The same interaction is consistently used to generate the 
mean-field, the excited states and the couplings. In this 
study, special attention is given to the issue of the double 
counting of the uncorrelated second-order diagram. The 
subtraction of this second-order term is shown not to lead 
to pathological behaviors when positive incident energy 
is considered, contrarily to what is expected in Ref. m- 
Moreover, the use of the finite-range Gogny interaction 
prevents from the necessity of ad-hoc momentum cut-off 
when second-order effects are considered. 


1.4 Project 

In the continuity of work presented in Ref. [46] , we wish to 
explore the possible connections provided by a microscopic 
optical potential using NSM. The framework is summa¬ 
rized in Fig. [3] In the first line of the diagram, effective 



Fig. 2. Connections. Dashed arrows emphasize connections 
discussed in the following. 


interaction is used within structure models to make spec¬ 
troscopic predictions. Feedback from experiment provides 


constraints on the interaction whenever a reliable struc¬ 
ture model is used. Effective interaction and structure cal¬ 
culation can then be used to define the microscopic po¬ 
tential through NSM. Once the corresponding scattering 
problem is solved, feedbacks are made possible from cross 
section data providing reaction constraints to the fit of 
the interaction. Microscopic potentials are nonlocal, com¬ 
plex and energy dependent. They can provide prescrip¬ 
tions for future phenomenological potentials, in particular 
concerning the shape of the nonlocality, the energy de¬ 
pendence. One can as well investigate the origin of the 
volume and the surface part of the potential in terms of 
target excitations. Reciprocally, when data are available, 
phenomenological potentials can help identifying contri¬ 
butions missing in microscopic potentials. Moreover when¬ 
ever phenomenological potentials obey dispersion relation, 
a connection is also made with spectroscopy. In Fig. [3] the 
main connections we wish to investigate in the following 
are highlighted. 

In the following, the NSM formalism for spherical- 
target nuclei and the integro-differential Schrodinger equa¬ 
tion are briefly exposed in Sec. 12.11 and Sec. 12.21 respec¬ 
tively. In Sec. 12.31 we emphasize the importance of the 
exact treatment of the intermediate HF propagator and 
more precisely the account of single-particle resonances. 
As a proof of concept in Sec. jS] we apply NSM to nucleon 
elastic scattering from "‘'’Ca. Some possible crosstalks be¬ 
tween phenomenological potentials and their microscopic 
counter part are discussed in Sec. |3| In Sec. [S] we show 
how phenomenological nonlocal potential can relate to the 
effective NN interaction through volume integrals. Finally 
in Sec. m we draw plans for the decade to come. In par¬ 
ticular, we mention the issue of spherical target nuclei ex¬ 
periencing pairing correlations and the one of deformed 
target nuclei. 

2 Method for spherical target 

2.1 NSM potential 

The NSM formalism is presented in detail in Ref. [13]. We 
briefly introduce the key points of the formalism. Equa¬ 
tions are presented omitting spin for simplicity. The po¬ 
tential, V, consists of two components, 

V = + AV. (1) 

This potential will be referred to as the NSM potential 
in the following. The former is a mean-field potential; the 
latter is a polarization potential built from target nucleus 
excitations. The HF potential, , is the major contri¬ 
bution to the real part of the optical potential. The polar¬ 
ization potential, AV, brings only a correction to the real 
part of V and entirely generates its imaginary contribu¬ 
tion. The HF potential is obtained in the Green’s function 
formalism neglecting two-body correlations m- It reads, 

I/^^(r, r’)= J driv{r,ri)p{ri)6{r — V) — v{r,V)p{r,r'), 

(2) 
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where v is the effective NN interaction and 

Pi^) = ( 3 ) 

i 

p(r,r’) = (4) 

i 

are the local and nonlocal densities with i running over oc¬ 
cupied states. is made of a local direct term and an 
exchange term which is nonlocal because of the finite range 
of Gogny interaction. It is energy independent. Rearrange¬ 
ment contributions stemming from the density-dependent 
term of the interaction are also taken into account. To¬ 
gether with Schrddinger equation, Eq. defines a self- 
consistent scheme as shown in Fig. [3] 


for once in and twice in When two-body cor¬ 

relations are neglected in Eq. one expects AV to re¬ 
duce to As a matter of fact shall be subtracted 
twice [42]. This formalism takes into account all the cor¬ 
relations explicitly. Although it is well suited for ab-initio 
developments, we wish to make the connection with effec¬ 
tive NN interactions. In practice, if one uses an effective 
interaction with a density-dependent term, such as Gogny 
or Skyrme forces, attention must be paid to correlations 
already accounted for in the interaction |40] . Indeed, in 
such a case, part of particle-particle correlations is already 
contained at the HE level. We thus use the same prescrip¬ 
tion as in Ref. [31], omitting the real part of while 
approximating the imaginary part of by Im [E^^^]. 
Then Eq. ([SJ reduces to 


NN interaction 



Schrodinger equation 
Fig. 3. Self-consistent Hartree-Fock. 


The polarization potential, AV in Eq. o, is built cou¬ 
pling the elastic channel to the intermediate excited states 
of the target nucleus. Those excited states are described 
within the RPA formalism. Both excited states and cou¬ 
plings are generated using Gogny interaction. Together 
with Schrodinger equation, Eq. © now defines a new self- 
consistent scheme, illustrated in Fig. |3| when considering 
both full-line and dashed-line arrows. In practice, as de- 


Z\E = Im 



+ yRPA 




( 6 ) 


Both ingredients of Eq. (jS]), and V^'^\ can be ex¬ 

pressed in terms of diagrams. In Fig. [5] wavy lines stand 
for the effective NN interaction and up (down) arrows 
stand for HE particle (hole) propagators. Subscript p {h) 
refers to the quantum numbers of the single-particle (hole) 
HE states used to build target excitations. The subscript A 
refers to the quantum numbers of the intermediate single¬ 
particle state of the scattered nucleon. Both discrete and 
continuum spectra of the intermediate single-particle state 
are accounted for. Label (a) refers to built with an 

unoccupied intermediate state. Label (a’) refers to 
built with an occupied intermediate state. Labels (b) and 
(b’) refer to the corresponding uncorrelated particle-hole 
contributions. Exchange diagrams are taken into account. 


NN interaction 



Schrodinger equation 


Fig. 4. Self-consistent RPA (full-line and dashed-lined ar¬ 
rows). Consistent RPA on top of a self-consistent HF (full-line 
arrow only). The HF propagator is dressed only once. 

scribed in Fig. UJ we first converge the HF scheme, as 
shown in Fig. (3] then we dress only once the HF prop¬ 
agator by coupling to the excitations of the target. This 
makes the scheme only consistent at that stage. Going into 
more details, the polarization contribution to the potential 

o p-| c 

AE = E^^-bE^^^-2E(2), (5) 

where and are contributions from particle- 

particle and particle-hole correlations, respectively. The 
uncorrelated particle-hole contribution, V^‘^\ is accounted 


(a) (b) 



Fig. 5. Diagrammatic contributions of (a) and E^^^ (b). 

Indices p, h and A refer to particle, hole and the intermediate 
state in the HF field, respectively. Wavy lines stand for the 
effective NN interaction. 
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For nucleons with incident energy the RPA potential, 
corresponding to Figs. [5^ and [5^’, reads, 





_^ 

E — ex-\- Eki 


ir{EN) 


E — Ex — Em + ir[EM) _ 

X n^{v)Q^{v% (7) 


where rii and Ei are occupation number and energy of the 
single-particle state (j)i in the HF field, respectively. Em 
and r{EM) represent the energy and the width of the 
excited state of the target, respectively. Additionally, 




(pA) 


where X and Y denote the usual RPA amplitudes and 


E^J\{r) 



ri(^*(ri)u(r,ri) 


1 


(l)xir)^j{ri), (9) 


where P is a particle-exchange operator and v is the same 
effective NN interaction as in Eq. The uncorrelated 
particle-hole contribution, corresponding to Figs. [SJd and 
133’, reads 


E(2)(r 


’ 2 

ij A 


n^{l - nj)nx 


E — Ex + Eij — ir(Eij) 


+ 


nj{l - ni){l - nx) 


E — Ex — Eij ir(Eij)_ 

X E,,x{r)E*^xir'), 


( 10 ) 


with Eij = Ei — Ej , the uncorrelated particle-hole energy. 


In practice, is determined in coordinate space to 
ensure the correct asymptotic behavior of single-particle 
states. This nonlocal potential is then used to build the 
single-particle intermediate state used to determine the 
polarization potential, AV. It is worth mentioning that 
the HF potential in coordinate space reproduces bound 
state energies obtained with the HF/DIS code on oscilla¬ 
tor basis which makes our result reliable. 

The description of target excitations is obtained solving 
RPA equations in a harmonic oscillator basis, including 
fifteen major shells [13] and using the Gogny DIS inter¬ 
action m- We account for RPA excited states with spin 
up to J = 8, including both parities, in order to achieve 
the convergence of cross section calculations. The first 
J'^ = 1“ excited state given by RPA, containing the spu¬ 
rious translational mode, is removed from the calculation. 
Moreover, in order to avoid spurious modes in the uncor¬ 
related particle-hole term, we approximate the = 1“ 
contribution in by half that of the = 1“ contri¬ 
bution in 

Even though RPA/DIS method provides a good over¬ 
all description of the spectroscopic properties of double- 
closed shell nuclei, still some contributions are left out. 


First, the projection on an oscillator basis discretizes the 
RPA continuum. As a consequence, the escape width is 
missing from the structure calculation. Second, couplings 
to two or more particle-hole states are excluded from the 
model space even though they may play a significant role. 

The impact of these couplings is a strength redistribution 
through a damping width as well as a shift in energy of ex¬ 
cited states. Third, the optical potential is, by definition, 
built to provide the energy-averaged iS-matrix. Hence, the 
rapid fluctuations, the potential exhibits at low energy 
due to compound-elastic contribution, shall be averaged 
before identifying the result of Eq. o with an optical po¬ 
tential [50] . 

In the present work, we simulate those three different 
widths assigning a single phenomenological width, T(Ajv), 
to each RPA state. E^Em) takes the value of 2, 5, 15 
and 50 MeV, for excitation energies of 20, 50, 100 and 
200 MeV, respectively. Those values have not been fitted 
in order to better reproduce cross section data. Longer 
term solutions are planned in order to provide more mi¬ 
croscopic prescriptions for those widths. The escape width 
can be obtained using continuum RPA |5T1H5] . We also 
plan to determine the damping width and the energy shift 
using the multiparticle-multihole configuration mixing method 

m 


2.2 Integro-differential Schrbdinger equation 


The integro-differential Schrodinger equation for bound 
states and scattering is solved without localization proce¬ 
dures. The radial Schrodinger equation reads. 


^ \f_ 

2fi dr'^ 


J.2 


fijir) 


+r 



'■,E)fi,ir'ydr' 


Efuir), (11) 


where fijir) = r<pij (r) is the partial wave for the projectile- 
target relative motion, E is the incident nucleon energy 
and vij(r,r'\ E) is defined from the multipole expansion 
of the nonlocal potential 


V(rcr,r’cr'; A) = ^ (r,/; F;)v/^.^(f'cr'), 

Ijm 

( 12 ) 

with 

yi]m {ra) = [Yi{r) (^Xi/ 2 {cr)]jm- (13) 

The potential, V, is complex and energy dependent for 
E > 0, and real and energy independent for A < 0. Dis¬ 
crete solutions are obtained by expressing Eq. m on a 
mesh in coordinate space and performing the correspond¬ 
ing matrix diagonalization [^. For positive energies, the 
scattering problem with the correct asymptotic conditions 
turns into a matrix inversion following J. Raynal’s method 
for scattering exposed in the DWBA code explanatory 
leaflet [Mj . 
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2.3 Resonances in the intermediate wave 

We now would like to emphasize the impact of an ex¬ 
act treatment of 4)\ on the second-order terms of the 
potential, and td(2) (ggg Eqs_ ([7]) through (IT?il) i: 

especially, the role of single-particle resonances already 
discussed by Rao et al. ES] within a phenomenological 
approach. In previous works, (j)\ has often been approxi¬ 
mated by a plane wave for neutron and a Coulomb wave 
for proton [40] or discretized [39] . In this work, we include 
both discrete and continuum spectra of 4>x determined in 
the HF field with the correct asymptotic solutions. Phase- 
shifts for neutron and proton scattering from ^'^Ca in the 
HF field are shown in Fig. [51 We observe single-particle 
resonances for several partial waves each time phaseshift 
increases rapidly through an odd multiple of 7r/2. Reso¬ 
nance energies are summarized in Table [TJ Those reso- 






Fig. 6. Phaseshift in the HF held for neutron (panel (a) and 
(b)) and proton (panel (c) and (d)) scattering from '^^Ga as a 
function of incident energy for the hrst ten partial waves. 


Table 1. HF single-particle resonance energies (in MeV) for 
neutron (n) and proton (p) scattering from "^^Ca. 



nances will result in fluctuations of the imaginary part of 
yRPA y( 2 ) Q and (fTOll 'l whenever the energy 

E—E]\f matches a resonance energy ex of the intermediate 
single-particle state. As a consequence, those resonances 
will strongly influence the corresponding cross section. As 



Fig. 7. Reaction cross section vs. proton incident energy with 
V = + Im[F^^^/2] coupling only to the hrst 1“ excited 

state in '*°Ca {Ex = 9.7 MeV). Comparison between (j)x treated 
as a Coulomb wave (dashed line) or as a single-particle state 
in the HF held (solid line). 


an illustration of resonance effects, we show, in Fig. [7] 
the reaction cross section for proton scattering with a po¬ 
tential including the HF potential as real part and the 
imaginary part of the RPA potential generated restricting 
couplings to the second = 1“ excited state of ^°Ca 
with Ex = 9.7 MeV excitation energy. This is done tak¬ 
ing T(Ajv) = 0 MeV in Eq. ([7]) in order to emphasize 
the effect of resonances. This result is compared with a 
calculation using a Coulomb wave as intermediate state. 
The exact treatment of the intermediate state leads to a 
global enhancement of the reaction cross section. Moreover 



Fig. 8. Reaction cross section vs. neutron incident energy. Neu¬ 
tron scattering from '‘°Ca with + Im (solid line) 

and + Im (dashed line) coupling to all available 

open channels. 
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we notice that coupling to only one excited state of the 
target already leads to four resonances between 10 and 20 
MeV. One expects to get a large number of resonant con¬ 
tributions once coupling to the thousand target excited 
states. As an example, we show in Fig. [5] the same cal¬ 
culation as in Fig. [7] but including all the open channels 
for a given incident energy. Once again as discussed in 
Sec. 12.11 one needs to average fluctuating contributions 
before comparing scattering observables with experiment. 
This shows the importance of a complete treatment of the 
intermediate wave. 


3 Microscopic potential and cross section 
data 

As a first application, NSM has been applied to neutron 
and proton scattering from ^°Ca using Gogny DIS inter¬ 
action. The corresponding differential cross sections for 
incident energies below 40 MeV are presented in Fig. |9l 
Compound-elastic corrections furnished by the Hauser- 
Feshbach formalism using Koning-Delaroche potential with 
TALYS are applied to cross sections obtained from NSM 
and Koning-Delaroche potential, respectively. It is mostly 
relevant below 10 MeV for neutron projectile while it gives 
a smaller contribution for proton. NSM results compare 
very well to experiment and those based on Koning-Delaroche 
potential up to about 30 MeV incident energy. References 
to data are given in Ref. [52] . Error bars are smaller than 
the size symbols. Beyond 30 MeV, backward-angle cross 
sections are overestimated. Discrepancies at 16.9 MeV (23.5 MeV) 
for neutron (proton) scattering are related to resonances 
in the intermediate single-particle state when not com¬ 
pletely averaged. A detailed treatment of the width might 
cure this issue. In Fig. (TU] we show calculated analyz¬ 
ing powers for neutron and proton scattering at several 
energies, in good agreement with measurements. More¬ 
over, agreement with data is comparable to that obtained 
from Koning-Delaroche potential. These results suggest 
that NSM potential retains the correct spin-orbit behav¬ 
ior. In Fig. 111! we show reaction cross section for proton 
scattering and total cross section for neutron scattering. 
Calculated reaction cross sections are in good agreement 
with experiments. For neutrons, however, we underesti¬ 
mate the total cross section below 10 MeV. Consider¬ 
ing that the differential elastic cross sections are well re¬ 
produced (see Fig. a), this underestimation suggests that 
part of absorption mechanisms is not accounted for, such 
as target-excited states beyond RPA, double-charge ex¬ 
change process or an intermediate deuteron formation. 

This microscopic potential makes the bridge between cross 
section data and effective interaction. It is worth mention¬ 
ing that we already get nice agreement with data without 
any adjustable parameter and using an effective interac¬ 
tion, the Gogny DIS interaction, originally tailored for 
structure purposes. This framework opens the way for new 
effective interactions based both on structure and reaction 
constraints. 



6 (deg.) 

c.m. ^ ® ^ 



Fig. 9. Differential cross sections for neutron (a) and proton 
(b) scattering from '*°Ca. Comparison between data (symbols), 
+ AV results (solid curves) and Koning-Delaroche poten¬ 
tial results (dashed curves). 
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(deg) 





0 (deg) 

c.m. ^ 

Fig. 10. Same as Fig. for analyzing powers. 

4 Phenomenological potential and 
microscopic potential 

We shall now highlight the possible crosstalk between mi¬ 
croscopic and phenomenological potentials depicted in Fig. 
A large variety of local potentials have been developed 
in order to describe reaction data. Mahaux and Sartor 
have then demonstrated the need for the potential to sat¬ 
isfy a dispersion relation, connecting its imaginary part 
to its real part, which provides a link with shell model 
[5S] . Along that line, local dispersive potentials have been 
developed [571158] . One issue in those local approaches is 
the spurious energy dependence of the potential coming 
from the use of a local ersatz to represent a nonlocal ob- 


Fig. 11. Reaction cross section for proton (a) and total cross 
section for neutron (b) scattering from ^°Ca. Gomparison be¬ 
tween data (symbols), -\- AV results (solid curve) and 

Koning-Delaroche potential (dashed curve). 


ject. This issue has been overcome building a dispersive 
potential with a nonlocal static real component (SU). A 
recent version of this dispersive potential is fully nonlocal 
both in its real and imaginary parts m- It is parametrized 
only for "‘^Ca but using all the structure and reaction data 
available for this nucleus. Another nonlocal dispersive po¬ 
tential is currently being developed for a broader range of 
nuclei m- It is interesting to compare such phenomeno¬ 
logical potentials with microscopic and ab-initio ones [3,81 
nsg. This connection can help identifying missing compo¬ 
nents in microscopic potentials. Reciprocally microscopic 
potential can provide some guidance for next-to-come po¬ 
tential parametrizations regarding for example the shape 
and the range of the nonlocality or the incident energy 
dependence. 

As an illustration, we now consider the case of neutron 
scattering from "‘^Ca at E = 9.91 MeV. In Fig. |T2l we 
compare the NSM potential with the nonlocal dispersive 
(NLD) potential from Ref. |5^. We focus on the multipole 
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expansion (see Eq. (IT^ l of the imaginary part of both 
nonlocal potentials and depict their diagonal contribu¬ 
tions. NSM and NLD potentials compare very well around 
r = 4.3 fm at the nucleus surface. In the volume region, 
NSM provides a much stronger contribution to the imag¬ 
inary potential than NLD potential. Nevertheless, NSM 
leads to a reasonable agreement with cross section data 
as shown in Fig. [HI At relatively low incident energy the 
projectile nucleon is only sensitive to the surface region. 
In other words, low-energy elastic observables do not con¬ 
strain the volume contribution of the potential. Further 
inelastic scattering calculations using both potentials may 
help disentangling this discrepancy as partial-wave func¬ 
tions play a major role in such calculations. In Fig ]131 we 


NLD - 

V«%AV - 



Fig. 12. Multipole expansion of the imaginary part of 
+ AV (solid curve) and NLD potential (dashed curve) for 
r = r', as a function of radius and partial waves for n-t-^°Ca 
at 9.91 MeV. 


compare the nonlocality of both imaginary components at 
the surface of the nucleus (r = 4.3 fm) as a function of the 
nonlocal parameter s = |r — r’|. We get a good agreement 
between the microscopic approach and the phenomenolog¬ 
ical one. Even though small emissive contributions appear 
in some of the multipoles, NSM validates the choice of a 
Gaussian nonlocality as originally proposed by Perey and 
Buck [521 used as well for NLD potential. NSM also re¬ 
produces the range of the nonlocality of the NLD potential 
which corresponds to a nonlocality parameter /3 = 0.94 fm 
at the surface of the nucleus. 

We have presented here results only for a given projec¬ 
tile at a given incident energy. A more exhaustive and 
hopefully conclusive study is in progress. Along the same 
line, it will be interesting to look at the nucleon asym¬ 
metry dependence of the NSM potential for example go¬ 
ing toward neutron-rich Ca isotopes and comparing with 
the dispersive potential obtained by Charity et al. for 

40,42,44,48,60,70Qg^ 



h 
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Fig. 13. Same as Fig. [12] for nonlocality at r = 4.3 fm. 


5 Phenomenological potential and effective 
NN interaction 

We now consider the connection between phenomenologi¬ 
cal potential and effective NN interaction through the mi¬ 
croscopic potential and the structure calculation as shown 
in Fig. [21 The link between microscopic potential and phe¬ 
nomenological one is done using volume integrals. Volume 
integrals are useful means of comparison between poten¬ 
tials as they are well constrained by scattering data. Here 
we focus on the real part of the potential. The volume 
integral for a given multipole (l,j) of the real part of the 
nonlocal potential is defined as, 

Jy{E) = J dr j dr'r'^Re[i^ij{r,r', E)], (14) 

where A is the nucleon number of the target. The HF po¬ 
tential is the leading contribution to the real part of NSM 
potential in Eq. O- In Fig. [TTJ we present the volume 
integral of the multipole expansion of the HF potential, 
Eq. ([2|), as a function of the partial wave. We compare it 
to the same quantity obtained from Perey-Buck (PB) non¬ 
local potential [52] • HF potential gives results similar to 
PB potential up to about the twelfth partial wave. Black 
segments denote the strongest partial-wave contributions 
accounting for 80% of the reaction cross section at the se¬ 
lected incident energies. Hence taking the PB potential as 
a reference, the HF potential has a reasonable behavior 
up to about 17 MeV incident energy. Beyond this energy 
HF saturates, following the trend of the Hartree poten¬ 
tial which is local and thus partial-wave independent. As 
a result increasing incident energy, HF yields a much too 
large volume integral reflected in an overestimate of the 
differential cross section at backward angles as shown in 
Fig. ini We present here results obtained with the DIS 
parametrization of the Gogny force, but we came to the 
same conclusions using the DIM parametrization [Sj. The 
behavior of the volume integral as a function of the partial 
wave is dictated by the shape and the range of the non¬ 
locality. PB potential is built with a Gaussian nonlocality 
whereas the HF potential is made of a local Hartree term 
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Fig. 14. Volume integral as a function of partial waves for neu¬ 
tron scattering from "^^Ga: HF potential (solid curve), Hartree 
potential (dash-dotted curve) and Perey-Buck potential (dot¬ 
ted curve). Horizontal segments denote the partial-wave inter¬ 
val to sum up 80% of the reaction cross section at selected 
incident energies. 


and a nonlocal Fock term as shown in Eq. Those two 
contributions can be related to the different terms of the 
effective NN interaction. Gogny interaction is built with 
two Gaussian ranges in its central part. The contributions 
to the Hartree potential of those two Gaussians as well 
as the contribution of the density-dependent term of the 
interaction are shown in panel (a) of Fig. fThl In panel (b), 
we show the contributions of the two central terms with 
finite range to the nonlocality of the first partial wave 
of the Fock term. The summation of those two Gaussian 
nonlocalities with opposite sign yields the ”W” shape of 
the total nonlocality. We then present in panel (c) the 
volume integral corresponding to the nonlocal Fock po¬ 
tential. The good behavior of the HF volume integral for 
low partial waves is related to the downward slope of the 
Fock volume integral and its change of sign around the 
seventh partial wave. Then for higher partial waves, the 
Fock volume integral converges to the zero limit. The cor¬ 
responding behavior of the Fock term as a function of the 
partial wave is depicted in panel (d). 

We plan to investigate to what extent the effective interac¬ 
tion could be improved in order to get a better behavior of 
the HF potential above 30 MeV. This issue emerges from 
the use of a local and energy-independent ersatz for the ef¬ 
fective NN interaction instead of the actual g-matrix with 
its full complexity. One way to get a fully nonlocal HF po¬ 
tential would be to use a nonlocal effective NN interaction 
as proposed for example by Tabakin |64] . 


6 Next decade 

The present work constitutes a promising step forward 
aimed to a model keeping at the same footing both re- 



0 5 10 15 20 25 30 


PW 



s (fm) 



s (fm) 


Fig. 15. Gontributions for n -|- ^'^Ga to: (a) to the Hartree local 
potential (V^): Total (solid line), first range of DIS (dashed 
line), second range of DIS (dash-dotted line) and density term 
(dotted line), (b) First partial wave of the nonlocal Fock term 
at r = r' = 4.3 fm: Total (solid line), first range of DIS (dashed 
line) and second range of DIS (dash-dotted line), (c) Volume 
integral of the Fock potential as a function of partial wave: Neg¬ 
ative slope (solid line), positive slope (dashed line), (d) Same as 
(c) for the Fock components nonlocality at r = r' = 4.3 fm. 


action and structure aspects of the many-nucleon sys¬ 
tem. Starting from an effective NN interaction, NSM ac¬ 
counts reasonably well for low-energy scattering data. We 
use consistently the Gogny DIS interaction, although this 
scheme can be applied to any interaction of similar nature. 
An important feature of the approach is the extraction of 
the imaginary part of the potential by means of intermedi¬ 
ate excitations of the target. The study has been restricted 
to closed-shell target but can be extended to account for 
pairing correlations as well as axial deformation. We now 
expose in more detail our plans for the forthcoming years. 


6.1 Spherical target nuclei 

In the short-term future, we wish to investigate the NSM 
scheme for spherical targets according to the following 
plan, 

. NSM will be applied to a broader range of target nuclei 
well described within RPA including "‘^Ga, ®°Zr, ^^^Sn 
and ^°®Pb. 

. The link between NSM potential and phenomenology 
can be carried on exploring energy, multipole and mass 
dependences of the potential. The NSM potential can 
also provide some trends for the shape and the range 
of the nonlocality. We plan to investigate the building 
of surface and volume contributions as a function of 
target excitations. 






















G. Blanchon et al.: Prospective study on microscopic potential with Gogny interaction 


11 


6.3 Deformed target nuclei 


Above about 50 MeV incident energy, a connection 
can be established between the NSM potential and the 
folding potential relying on g-matrix and thus with the 
bare NN interaction. This can be a fertile ground for 
new effective interactions. 

At low incident energy, NSM provides a volume part 
of the imaginary potential larger than phenomenology. 

The small volume contribution in phenomenological 
potential is often justihed by the fact that the pro¬ 
jectile does not have sufficient energy to knock out a 
target nucleon. This discrepancy is possibly due to the 
fact that at low energy the volume part of the imagi¬ 
nary potential is not well constrained because the pro¬ 
jectile does not explore the interior of the potential. We 
plan to use the NSM potential in inelastic scattering 
calculations in order to disentangle this issue. 

In its present version, NSM requires a phenomeno¬ 
logical width. This width has several microscopic ori¬ 
gins, as discussed in Sec. 12.11 A microscopic account 
of those widths is planed using continuum RPA m 
and multihole-multiparticle conhguration mixing [52) 
for the escape and the damping widths, respectively. 

The NSM potential will be used to provide transmis¬ 
sion coefficients for compound-elastic calculations. More- under Grant No 1120396. 
over, we plan to develop a compound-nucleus formal¬ 
ism based only on NSM. Indeed, NSM gives access to 
the fluctuating contribution of the S-matrix and as a 
consequence to the compound-elastic contribution. 

The study of the volume integral of the real part of the 
potential has exhibited the possible crosstalk between 
phenomenological potentials and effective NN interac¬ 
tions. In particular, the interaction is not well suited 
for the description of partial waves with more than 
about £ = 7. Those reaction constraints will be used 
for new parametrizations of the interaction. Moreover, 
a nonlocal version of the effective NN interaction could 
tackle the issue of the saturation of the HF volume in¬ 
tegral to the Hartree one. 


In the midterm future, we plan to deal with axially-deformed 
targets according to the following plan, 

. This will require the development of an axially-deformed 
HFB potential in coordinate space. The correspond¬ 
ing mean-held and the pairing held have already been 
studied in coordinate space in a previous work on Cooper’s 
pairs [55] . 

. QRPA [2 will be used to generate excited states in the 
intrinsic frame of the target. 

• A projection an good angular momentum, using the ro¬ 
tational approximation [70] , will provide the monopole 
and different coupling potentials to model nucleon elas¬ 
tic scattering from axially deformed target. 

. The problem of solving coupled equations with nonlo¬ 
cal potentials will have to be addressed. 
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